Transverse momentum dependence of 77 meson suppression 
in Au+Au collisions at ^/snn = 200 GeV 

A. Adare, 12 S. Afanasiev, 27 C. Aidala, 40 N.N. Ajitanand, 57 Y. Akiba, 51 ' 52 H. Al-Bataineh, 46 J. Alexander, 57 
K. Aoki, 33 ' 51 L. Aphecetche, 59 Y. Aramaki, 11 J. Asai, 51 E.T. Atomssa, 34 R. Averbeck, 58 T.C. Awes, 47 B. Azmoun, 6 
V. Babintsev, 23 M. Bai, 5 G. Baksay, 19 L. Baksay, 19 A. Baldisseri, 15 K.N. Barish, 7 P.D. Barnes, 36 B. Bassalleck, 45 
A.T. Basye, 1 S. Bathe, 7 S. Batsouli, 47 V. Baublis, 50 C. Baumann, 41 A. Bazilevsky, 6 S. Belikov, 6 'B R- Belmont, 63 
R. Bennett, 58 A. Berdnikov, 54 Y. Berdnikov, 54 A. A. Bickley, 12 J.G. Boissevain, 36 J.S. Bok, 66 H. Borel, 15 K. Boyle, 58 

M.L. Brooks, 36 H. Buesching, 6 V. Bumazhnov, 23 G. Bunce, 6 ' 52 S. Butsyk, 36 CM. Camacho, 36 S. Campbell, 58 
B.S. Chang, 66 W.C. Chang, 2 J.-L. Charvet, 15 C.-H. Chen, 58 S. Chernichenko, 23 CY. Chi, 13 M. Chiu, 6 ' 24 I.J. Choi, 66 
R.K. Choudhury, 4 P. Christiansen, 38 T. Chujo, 62 P. Chung, 57 A. Churyn, 23 O. Chvala, 7 V. Cianciolo, 47 Z. Citron, 58 
B.A. Cole, 13 M. Connors, 58 P. Constantin, 36 M. Csanad, 17 T. Csorgo, 30 T. Dahms, 58 S. Dairaku, 33 ' 51 I. Danchev, 63 
K. Das, 20 A. Datta, 40 G. David, 6 A. Denisov, 23 D. d'Enterria, 34 A. Deshpande, 52 ' 58 E.J. Desmond, 6 O. Dietzsch, 55 

A. Dion, 58 M. Donadelli, 55 O. Drapier, 34 A. Drees, 58 K.A. Drees, 5 A.K. Dubey, 65 J.M. Durham, 58 A. Durum, 23 
D. Dutta, 4 V. Dzhordzhadze, 7 S. Edwards, 20 Y.V. Efremenko, 47 F. Ellinghaus, 12 T. Engelmore, 13 A. Enokizono, 35 

H. En'yo, 51 ' 52 S. Esumi, 62 K.O. Eyser, 7 B. Fadem, 42 D.E. Fields, 45 ' 52 M. Finger, Jr., 8 M. Finger, 8 F. Fleuret, 34 
S.L. Fokin, 32 Z. Fraenkel, 65 'B J.E. Frantz, 58 A. Franz, 6 A.D. Frawley, 20 K. Fujiwara, 51 Y. Fukao, 33 ' 51 T. Fusayasu, 44 
I. Garishvili, 60 A. Glenn, 12 H. Gong, 58 M. Gonin, 34 J. Gosset, 15 Y. Goto, 51 ' 52 R. Granier de Cassagnac, 34 
N. Grau, 13 S.V. Greene, 63 M. Crosse Perdekamp, 24 ' 52 T. Gunji, 11 H.-A. Gustafsson, 38 'B A. Hadj Henni, 59 
J.S. Haggerty, 6 K.I. Hahn, 18 H. Hamagaki, 11 J. Hamblen, 60 J. Hanks, 13 R. Han, 49 E.P. Hartouni, 35 K. Haruna, 22 

E. Haslum, 38 R. Hayano, 11 M. Heffner, 35 S. Hegyi, 30 T.K. Hemmick, 58 T. Hester, 7 X. He, 21 J.C. Hill, 26 
M. Hohlmann, 19 W. Holzmann, 13 ' 57 K. Homma, 22 B. Hong, 31 T. Horaguchi, 11 ' 22 ' 51 ' 61 D. Hornback, 60 S. Huang, 63 

T. Ichihara, 51 ' 52 R. Ichimiya, 51 J. Ide, 42 H. Iinuma, 33 ' 51 Y. Ikeda, 62 K. Imai, 33 * 51 J. Imrek, 16 M. Inaba, 62 
D. Isenhower, 1 M. Ishihara, 51 T. Isobe, 11 M. Issah, 57 ' 63 A. Isupov, 27 D. Ivanischev, 50 B.V. Jacak, 58 'Q J. Jia, 6 ' 13 ' 57 

J. Jin, 13 B.M. Johnson, 6 K.S. Joo, 43 D. Jouan, 48 D.S. Jumper, 1 F. Kajihara, 11 S. Kametani, 51 N. Kamihara, 52 
J. Kamin, 58 J.H. Kang, 66 J. Kapustinsky, 36 D. Kawall, 40 ' 52 M. Kawashima, 53 ' 51 A.V. Kazantsev, 32 T. Kempel, 26 
A. Khanzadeev, 50 K.M. Kijima, 22 J. Kikuchi, 64 B.I. Kim, 31 D.H. Kim, 43 D.J. Kim, 28 ' 66 E.J. Kim, 9 E. Kim, 56 

S.H. Kim, 66 Y.J. Kim, 24 E. Kinney, 12 K. Kiriluk, 12 A. Kiss, 17 E. Kistenev, 6 J. Klay, 35 C. Klein-Boesing, 41 
L. Kochenda, 50 B. Komkov, 50 M. Konno, 62 J. Koster, 24 D. Kotchetkov, 45 A. Kozlov, 65 A. Krai, 14 A. Kravitz, 13 
G.J. Kunde, 36 K. Kurita, 53 ' 51 M. Kurosawa, 51 M.J. Kweon, 31 Y. K won, 60 ' 66 G.S. Kyle, 46 R. Lacey, 57 Y.S. Lai, 13 
J.G. Lajoie, 26 D. Layton, 24 A. Lebedev, 26 D.M. Lee, 36 J. Lee, 18 K.B. Lee, 31 K. Lee, 56 K.S. Lee, 31 T. Lee, 56 

M.J. Leitch, 36 M.A.L. Leite, 55 E. Leitner, 63 B. Lenzi, 55 P. Liebing, 52 L.A. Linden Levy, 12 T. Liska, 14 
A. Litvinenko, 27 H. Liu, 36 ' 46 M.X. Liu, 36 X. Li, 10 B. Love, 63 R. Luechtenborg, 41 D. Lynch, 6 C.F. Maguire, 63 
Y.I. Makdisi, 5 A. Malakhov, 27 M.D. Malik, 45 V.I. Manko, 32 E. Mannel, 13 Y. Mao, 49 - 51 L. Masek, 8 ' 25 H. Masui, 62 
F. Matathias, 13 M. McCumber, 58 P.L. McGaughey, 36 N. Means, 58 B. Meredith, 24 Y. Miake, 62 A.C. Mignerey, 39 

P MiVoS 8,25 v MiVi 62 A T\/m„, r 6 m lu^cVir-a 3 tt iu^KoU 6 A XT A/IWhanf-ir 4 v TU^r-;™ 11 A l\/Tnr™alD 7 



F. Matathias, 13 M. McCumber, 58 P.L. McGaughey, 36 N. Means, 58 B. Meredith, 24 Y. Miake, 62 A.C. Mignerey, 1 
P. Mikes, 8 ' 25 K. Miki, 62 A. Milov, 6 M. Mishra, 3 J.T. Mitchell, 6 A.K. Mohanty, 4 Y. Morino, 11 A. Morreale, 7 

D.P. Morrison, 6 T.V. Moukhanova, 32 D. Mukhopadhyay, 63 J. Murata, 53 ' 51 S. Nagamiya, 29 J.L. Nagle, 12 
M. Naglis, 65 M.I. Nagy, 17 I. Nakagawa, 51 ' 52 Y. Nakamiya, 22 T. Nakamura, 22 ' 29 K. Nakano, 51 ' 61 J. Newby, 35 
M. Nguyen, 58 T. Niita, 62 R. Nouicer, 6 A.S. Nyanin, 32 E. O'Brien, 6 S.X. Oda, 11 CA. Ogilvie, 26 K. Okada, 52 
M. Oka, 62 Y. Onuki, 51 A. Oskarsson, 38 M. Ouchida, 22 K. Ozawa, 11 R. Pak, 6 A.P.T. Palounek, 36 V. Pantuev, 58 
V. Papavassiliou, 46 I.H. Park, 18 J. Park, 56 S.K. Park, 31 W.J. Park, 31 S.F. Pate, 46 H. Pei, 26 J.-C. Peng, 24 

H. Pereira, 15 V. Peresedov, 27 DYu. Peressounko, 32 C. Pinkenburg, 6 R.P. Pisani, 6 M. Proissl, 58 M.L. Purschke, 6 
A.K. Purwar, 36 H. Qu, 21 J. Rak, 28 ' 45 A. Rakotozafmdrabe, 34 I. Ravinovich, 65 K.F. Read, 47 ' 60 S. Rembeczki, 19 
K. Reygers, 41 V. Riabov, 50 Y. Riabov, 50 E. Richardson, 39 D. Roach, 63 G. Roche, 37 S.D. Rolnick, 7 M. Rosati, 26 
CA. Rosen, 12 S.S.E. Rosendahl, 38 P. Rosnet, 37 P. Rukoyatkin, 27 P. Ruzicka, 25 V.L. Rykov, 51 B. Sahlmueller, 41 

N. Saito, 29 ' 33 ' 51 - 52 T. Sakaguchi, 6 S. Sakai, 62 K. Sakashita, 51 ' 61 V. Samsonov, 50 S. Sano, 11 ' 64 T. Sato, 62 
S. Sawada, 29 K. Sedgwick, 7 J. Seele, 12 R. Seidl, 24 AYu. Semenov, 26 V. Semenov, 23 R. Seto, 7 D. Sharma, 65 

I. Shein, 23 T.-A. Shibata, 51 ' 61 K. Shigaki, 22 M. Shimomura, 62 K. Shoji, 33 ' 51 P. Shukla, 4 A. Sickles, 6 C.L. Silva, 55 
D. Silvermyr, 47 C. Silvestre, 15 K.S. Sim, 31 B.K. Singh, 3 CP. Singh, 3 V. Singh, 3 M. Slunecka, 8 A. Soldatov, 23 
RA. Soltz, 35 W.E. Sondheim, 36 S.P. Sorensen, 60 I.V. Sourikova, 6 NA. Sparks, 1 F. Staley, 15 P.W. Stankus 47 

E. Stenlund, 38 M. Stepanov, 46 A. Ster, 30 S.P. Stoll, 6 T. Sugitate, 22 C. Suire, 48 A. Sukhanov, 6 J. Sziklai, 30 
E.M. Takagui, 55 A. Taketani, 51 ' 52 R. Tanabe, 62 Y. Tanaka, 44 K. Tanida, 33 ' 51 ' 52 ' 56 M.J. Tanncnbaum, 6 S. Tarafdar, 3 
A. Taranenko, 57 P. Tarjan, 16 H. Thcmann, 58 T.L. Thomas, 45 M. Togawa, 33 ' 51 A. Toia, 58 L. Tomasek, 25 Y. Tomita, 62 



2 



H. Torii, 22 ' 51 R.S. Towell, 1 V-N. Tram, 34 I. Tserruya, 65 Y. Tsuchimoto, 22 C. Vale, 6 - 26 H. Valle, 63 H.W. van Hecke, 36 
E. Vazqucz-Zambrano, 13 A. Vcicht, 24 J. Velkovska, 63 R. Vertesi, 16 - 30 A.A. Vinogradov, 32 M. Virius, 14 V. Vrba, 25 
E. Vznuzdaev, 50 X.R. Wang, 46 D. Watanabe, 22 K. Watanabe, 62 Y. Watanabe, 51 ' 52 F. Wei, 26 R. Wei, 57 
J. Wessels, 41 S.N. White, 6 D. Winter, 13 J.R Wood, 1 C.L. Woody, 6 R.M. Wright, 1 M. Wysocki, 12 W. Xie, 52 
Y.L. Yamaguchi, 11 ' 64 K. Yamaura, 22 R. Yang, 24 A. Yanovich, 23 J. Ying, 21 S. Yokkaichi, 51 - 52 G.R. Young, 47 
I. Younus 45 Z. You 49 I.E. Yushmanov, 32 W.A. Zajc, 13 O. Zaudtke, 41 C. Zhang, 47 S. Zhou, 10 and L. Zolin 27 

(PHENIX Collaboration) 

'Abilene Christian University, Abilene, Texas 79699, USA 
2 Institute of Physics, Academia Sinica, Taipei 11529, Taiwan 
'Department of Physics, Banaras Hindu University, Varanasi 221005, India 
4 Bhabha Atomic Research Centre, Bombay 400 085, India 
5 Collider-Accelerator Department, Brookhaven National Laboratory, Upton, New York 11973-5000, USA 
6 Physics Department, Brookhaven National Laboratory, Upton, New York 11973-5000, USA 
7 University of California - Riverside, Riverside, California 92521, USA 
8 Charles University, Ovocny trh 5, Praha 1, 116 36, Prague, Czech Republic 
9 Chonbuk National University, Jeonju 561-756, Korea 
10 China Institute of Atomic Energy (CIAE), Beijing, People's Republic of China 
11 Center for Nuclear Study, Graduate School of Science, University of Tokyo, 7-3-1 Hongo, Bunkyo, Tokyo 113-0033, Japan 

12 University of Colorado, Boulder, Colorado 80309, USA 
1:1 Columbia University, New York, New York 10027 and Nevis Laboratories, Irvington, New York 10533, USA 
14 Czech Technical University, Zikova 4, 166 36 Prague 6, Czech Republic 
15 Dapnia, CEA Saclay, F-91191, Gif-sur-Yvette, France 
16 Debrecen University, H-4010 Debrecen, Egyetem ter 1, Hungary 
11 ELTE, Eotv'ds Lordnd University, H - 1117 Budapest, Pdzmdny P. s. 1/A, Hungary 
18 Ewha Womans University, Seoul 120-750, Korea 
19 Florida Institute of Technology, Melbourne, Florida 32901, USA 
20 Florida State University, Tallahassee, Florida 32306, USA 
21 Georgia State University, Atlanta, Georgia 30303, USA 
22 Hiroshima University, Kagamiyama, Higashi- Hiroshima 739-8526, Japan 
23 IHEP Protvino, State Research Center of Russian Federation, Institute for High Energy Physics, Protvino, 142281, Russia 
24 University of Illinois at Urbana- Champaign, Urbana, Illinois 61801, USA 
2j Institute of Physics, Academy of Sciences of the Czech Republic, Na Slovance 2, 182 21 Prague 8, Czech Republic 

26 Iowa State University, Ames, Iowa 50011, USA 
27 Joint Institute for Nuclear Research, 141980 Dubna, Moscow Region, Russia 
28 Helsinki Institute of Physics and University of Jyvdskyld, P.O.Box 35, FI-40014 Jyvaskyla, Finland 
29 KEK, High Energy Accelerator Research Organization, Tsukuba, Ibaraki 305-0801, Japan 
:l " KFKI Research Institute for Particle and Nuclear Physics of the Hungarian Academy 
of Sciences (MTA KFKI RMKI), H-1525 Budapest 114, POBox 49, Budapest, Hungary 
31 Korea University, Seoul 136-701, Korea 
32 Russian Research Center "Kurchatov Institute", Moscow, Russia 
ss Kyoto University, Kyoto 606-8502, Japan 
34 Laboratoire Leprince-Ringuet, Ecole Polytechnique, CNRS-IN2P3, Route de Saclay, F-91128, Palaiseau, France 
35 Lawrence Livermore National Laboratory, Lwermore, California 94550, USA 
36 Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 
37 LPC, Universite Blaise Pascal, CNRS-IN2P3, Clermont-Fd, 63177 Aubiere Cedex, France 
38 Department of Physics, Lund University, Box 118, SE-221 00 Lund, Sweden 
39 University of Maryland, College Park, Maryland 20742, USA 
40 Department of Physics, University of Massachusetts, Amherst, Massachusetts 01003-9337, USA 
41 Institut fur Kernphysik, University of Muenster, D-48149 Muenster, Germany 
42 Muhlenberg College, Allentown, Pennsylvania 18104-5586, USA 
43 Myongji University, Yongin, Kyonggido 449-728, Korea 
44 Nagasaki Institute of Applied Science, Nagasaki-shi, Nagasaki 851-0193, Japan 
45 University of New Mexico, Albuquerque, New Mexico 87131, USA 
46 New Mexico State University, Las Cruces, New Mexico 88003, USA 
47 Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 
48 IPN-Orsay, Universite Paris Sud, CNRS-IN2P3, BP1, F-91406, Orsay, France 
4 Peking University, Beijing, People's Republic of China 
S0 PNPI, Petersburg Nuclear Physics Institute, Gatchina, Leningrad region, 188300, Russia 
51 RIKEN Nishma Center for Accelerator-Based Science, Wako, Saitama 351-0198, JAPAN 
S2 RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, New York 11973-5000, USA 
53 Physics Department, Rikkyo University, 3-34-1 Nishi-Ikebukuro, Toshima, Tokyo 171-8501, Japan 
54 Saint Petersburg State Polytechnic University, St. Petersburg, Russia 



3 



55 Universidade de Sao Paulo, Institute de Fisica, Caixa Postal 66318, Sao Paulo CEP05315-970, Brazil 
56 Seoul National University, Seoul 151-742, Korea 
5 Chemistry Department, Stony Brook University, Stony Brook, SUNY, New York 11794-3400, USA 
58 Department of Physics and Astronomy, Stony Brook University, SUNY, Stony Brook, New York 11794, USA 
59 SUBATECH (Ecole des Mines de Nantes, CNRS-IN2P3, Universite de Nantes) BP 20722 - 44307, Nantes, France 

60 University of Tennessee, Knoxville, Tennessee 37996, USA 
61 Department of Physics, Tokyo Institute of Technology, Oh-okayama, Meguro, Tokyo 152-8551, Japan 
62 Institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305, Japan 
es Vanderbilt University, Nashville, Tennessee 37235, USA 
6J * Waseda University, Advanced Research Institute for Science and 
Engineering, 17 Kikui-cho, Shinjuku-ku, Tokyo 162-0044, Japan 
65 Weizmann Institute, Rehovot 76100, Israel 
a6 Yonsei University, IPAP, Seoul 120-749, Korea 
(Dated: May 27, 2010) 

New measurements by the PHENIX experiment at RHIC for rj production at midrapidity as a 
function of transverse momentum (pr) and collision centrality in -y/ijvjv = 200 GeV Au+Au and 
p + p collisions are presented. They indicate nuclear modification factors (Raa) which are similar 
both in magnitude and trend to those found in earlier it measurements. Linear fits to -Raa as 
a function of pr in 5-20 GeV /c show that the slope is consistent with zero within two standard 
deviations at all centralities although a slow rise cannot be excluded. Having different statistical and 
systematic uncertainties, the ir° and r] measurements are complementary at high pr; thus, along 
with the extended pr range of these data they can provide additional constraints for theoretical 
modeling and the extraction of transport properties. 

PACS numbers: 25.75.Dw, 13.85.Qk, 13.20.Fc, 13.20.He 



Suppression of high pr hadron production in Au+Au 
collisions at RHIC [l|, [2| and its absence in <i+Au col- 
lisions [|[ provided the first direct evidence that an ex- 
tremely dense medium is formed in heavy ion collisions 
at RHIC energies. This suppression relative to the yield 
expected from the convolution of independent nucleon- 
nucleon scatterings, measured by the nuclear modifica- 
tion factor Raa, is now confirmed up to 20 GeV/c 
with identified tt° and attributed to the energy loss of 
the hard scattered partons in the dense medium. Sev- 
eral models with very different assumptions describe the 
magnitude of the observed tt° suppression, but predict 
slightly different evolution with increasing px- Calcu- 
lations based on perturbative QCD (pQCD) and static 
plasma predict that the fractional parton energy loss de- 
creases with pt like log{px) / px leading to a slow rise of 
the Raa with px (for a recent review see [U). In contrast, 
some AdS/CFT calculations find that the fractional en- 
ergy loss is proportional to px- Therefore,i?AA decreases 
with increasing transverse momentum [5j-l8| . The univer- 
sal upper bound model [t| predicts that Raa remains 
almost independent of the energy of the original gluon 
or quark. Other effects (modified nuclear parton distri- 
bution functions, Cronin-effect, modified fragmentation 
functions, the quark/gluon ratio) at given xx (2px/y/s) 
can also dependence of Raa > and it is clear change the px 
dependence of Raa, and it is clear that a precise mea- 
surement of the evolution of Raa with px would help 
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in confirming or rejecting classes of theories and putting 
tight constraints on the free parameters of the remaining 
ones. The first rigorous attempt to confront the observed 
7T° suppression with various pQCD-based parton energy 
loss calculations and to put quantitative constraints on 
the transport properties of the medium was made in [Toj 
using PHENIX 7r° data. One intriguing result was that a 
linear fit with a slope consistent with zero described the 
evolution of Raa with px slightly better than any of the 
pQCD models predicting a slow rise. However, the large 
statistical and systematic uncertainties of the high pt 7r° 
points prevented a clear distinction between constant or 
slowly rising Raa- 

One reason the ir° data Q allow such ambiguous in- 
terpretations is that the experimental uncertainties rise 
rapidly as we move to higher pt (>12-14 GeV/c), due 
to "shower merging," as explained below. In the case of 
the rj this problem is absent for px up to 50 GeV/c, sig- 
nificantly beyond the pt range expected to be accessible 
at RHIC. While the yield of the actually reconstructed rj 
mesons is smaller except at the highest px, the improve- 
ment in systematic uncertainties can help provide better 
constraints in comparisons to theory at high px and thus 
complement the ir° results. Of course some caution in in- 
terpreting the results is warranted: while both ir° and rj 
consist of light quarks, rj does have a hidden strangeness 
(ss) content so it is not a priori obvious that the ir° and 
r\ results are interchangeable. Earlier measurements [Til ] 
have shown that at least up to 12 GeV/c, the ir° and 
t] nuclear modification factors in Au+Au agree within 
uncertainties and the rj/ir ratio is constant for pt>4 
GeV/c in p+p ll]. Using recent, more precise measure- 
ments in PHENIX, we will re-examine whether ir° and 
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TABLE I: Typical systematic uncertainties on r\ spectra and 
Raa- See text for explanation of error types. 
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FIG. 1: (Color online) 77 invariant mass distribution for two 
different pair px bins (minimum bias data). Top: 9<pt<10 
GeV/c, the combinatorial background has been subtracted 
by using mixed events. Note the large difference between n° 
and f] raw yields. Insert: the r\ region magnified. Bottom: 
16<pr<18 GeV/c region, where mixed event subtraction is 
no longer necessary. Also, here a cut on the 7-pair energy 
asymmetry, a < 0.6 has been applied, which greatly improves 
the signal/background ratio at the rj peak but cuts into the 
lower part of the ty° peak due to cluster merging. 



77 production at midrapidity is indeed similar and study 
the asymptotic behavior of -Raa- 

This analysis used 3.25B minimum bias (MB) ^/snn = 
200 GeV Au+Au events, corresponding to 0.511 nb _1 
recorded in 2007 as well as 429M minimum bias (18.7 
nb- 1 ) and 2.06B triggered (6.90 pb" 1 ) y/s = 200 GeV 
p+p events recorded in 2006 in the PHENIX experiment 
at RHIC. Both the Au+Au and p+p data sets were ana- 
lyzed using the same analysis chain and cuts; thus, some 
of the systematic uncertainties cancel when we calculate 
the nuclear modification factor -Raa for Au+Au. 

Collision centrality in Au+Au has been established by 
the beam-beam counters [H[ (BBC, 3.0 < \t)\ < 3.9). A 
Glauber-model Monte Carlo along with a simulation of 
the BBC response was used to estimate the average num- 
ber of participating nucleons (iV par t) and binary nucleon- 
nucleon collisions (iVcoii) for each centrality bin 12 1. 

The 77 mesons were measured via their 77 — ¥ 77 decay 
channel. The photons were reconstructed in the lead- 
scintillator (PbSc) sectors of the PHENIX Electromag- 
netic Calorimeter (EMCal) 1. lj covering 3/8 of the full 
azimuth and —0.35 < 77 < 0.35 in pseudorapidity, and the 
77 yield was extracted from two-photon invariant mass dis- 
tributions. This analysis is similar to the one described 
in (Til [l5| . There are three important differences. In 
the case of ir° starting around pr=12 GeV/c the min- 
imum opening angle of the two decay photons is small 
enough for the photon showers to merge and become in- 
distinguishable. As pr increases, this effect leads to an 
increasing loss of observed 7r°, resulting in large correc- 
tions and corresponding systematic uncertainties (which 
are in fact the dominant systematic uncertainties at high 
Pt)- Since the mass of the 77 is about four times larger 
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than the tt° , this is not a problem for the 77 measurement 
up to p<r~50 GeV/c. On the other hand the observable 
?7 rates are much lower at low and medium px, as seen 
in the invariant mass distributions in Fig. [TJ because of 
the smaller branching ratio into two photons (39%) and 
the small ?7/7r «0.5 production ratio. The raw yields be- 
come comparable only around 20 GeV/c. Finally, in 
the 77 analysis we applied an a < 0.6 photon pair energy 
asymmetry cut (as opposed to a < 0.8 for 7r°) in order 
to improve the signal/background ratio in the r\ region. 



TABLE II: Parameters of the power-law fits A/p T for Au+Au 
and p+p . The errors used for fit are the statistical and 
pr-uncorrelated (Type A) systematic uncertainties added in 
quadrature. The pr range of the fits is 5-22 GeV/c. 



System/Cent. 
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n 


X 2 /NDF 


Au+Au 0-5 % 


27.2+11.9 


7.90+0.22 


3.1/7 


Au+Au 0-10 % 


17.6+5.5 


7.77+0.15 


10.6/8 


Au+Au 10-20 % 


19.1+5.9 


7.89+0.16 


10.2/9 


Au+Au 0-20 % 


18.5+4.3 


7.84+0.12 


10.5/7 


Au+Au 20-40 % 


17.3+4.2 


8.01+0.12 


17.2/8 


Au+Au 40-60 % 


9.53+2.65 


8.05+0.15 


5.5/8 


Au+Au 20-60 % 


14.5+2.5 


8.07+0.08 


11.2/9 


Au+Au 60-92 % 


1.13+0.40 


7.78+0.18 


2.98/6 


Au+Au MinBias 


10.4+1.4 


8.04+0.08 


9.41/9 


p+p 


8.84+0.99 


8.21+0.05 


8.33/9 



The raw 77 yield is always counted by integrating the 
histogram bin content in the 77 mass window (typically 
±30 MeV/c 2 ), but the way we treat the underlying com- 
binatorial background varies as a function of pr- In 
Au+Au up to 10 GeV/c, mixed event subtraction is 
used. The 77 region is then fitted with a polynomial and 
Gaussian (see insert in Fig. [T]) to estimate the residual 
background. When the signal/background ratio reaches 
1.0, already in the 7-10 GeV/c range, depending on 
centrality, mixed event subtraction is no longer needed; 
a polynomial and Gaussian fit is used on the original in- 
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FIG. 2: (Color online) Cross section of p + p — > r\ + X from 
the 2006 p+p data set (solid circles) and rj invariant yield in 
Au+Au collisions of various centralities (open symbols) and 
minimum bias (solid squares) from the 2007 data set. p+p 
is shown at the true pr value, all other spectra are shifted 
alternately by ±0.1 GeV/c for better visibility of the error 
bars and upper limits. 



variant mass distribution to estimate the background. At 
even higher pr (12-16 GeV/c) we estimate the residual 
background under the peak simply from the average bin 
content of the sidebands (the regions above and below 
the peak). 

Systematic uncertainties are classified into three types: 
Type A is pr-uncorrelated ( "point- by-point" ) and for the 
purposes of fitting and plotting, is added in quadrature 
to the statistical errors. Type C is the overall normaliza- 
tion uncertainty allowing all points to move by the same 
fraction up or down. Type B is all other ^-correlated 
uncertainties (including the cases where the shape of the 
correlation function is not known). Table fl] lists typi- 
cal uncertainties on the spectra and Raa- "Conversion 
(HBD)" stands for loss due to photon conversion in the 
Hadron Blind Detector, which was present in one of the 
two central arms during the 2007 (Au+Au) data taking. 
"ERT norm." stands for the normalization uncertainty of 
the EMCal-RICH Trigger, selecting high pr photons and 
electrons. "Acceptance variations" are small day-by-day 
changes of dead areas in the detector and thus are inde- 
pendent for the p+p and Au+Au runs. The systematic 
uncertainties on raw yield, photon PID and conversion 
(other) are common in p+p and Au+Au, and hence were 
partially cancelled out in the Raa calculation. 
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FIG. 3: (Color online) Nuclear modification factor for 77 
at various centralities, calculated using the measured p+p 
points. Dark (green) band around 1 indicates the absolute 
normalization error from p+p, light (grey) band is the (cen- 
trality dependent) absolute normalization error from Au+Au. 
Error bars include statistical and pr-uncorrelated systematic 
errors. Also shown: linear fits to the data with la error bands. 



Cross sections for p+p — > r\ + X and invariant yield 
of inclusive rj production in Au+Au collisions for dif- 
ferent centralities are shown in Fig. [2] They cover the 
5 < pt < 22 GeV/c range and five orders of magnitude 
in cross section (invariant yield). The overall normal- 
ization uncertainties (Type C) are 13% for p+p and 5% 
for Au+Au. Parameters of simple power-law fits (A/p T ) 
to various, partially overlapping centrality selections, in- 
cluding ones not shown in Fig. [2] are given in Table [TTJ 
Fits include all available points in the 5 <pr< 22 GeV/c 
range but exclude upper limits. Only statistical and pt- 
uncorrelated uncertainties were used in the fits. Note 
that for 7T° in Au+Au collisions the power n was con- 
sistent within uncertainties at all centralities ranging 
from 8.00+0.12 in 0-5% to 8.06+0.08 in 80-92%, and 
for 7T° in p+p the power n was 8.22+0.09. In this mea- 
surement we find that for 77 production p + p — > rj + X 
the power n is the same as it was for ir°. The powers 
obtained for rj in Au+Au are also consistent with those 
from 7T° within two standard deviations. 
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FIG. 4: (Color online) Nuclear modification factor -Raa for 
7T° (open squares, points shifted for clarity, data from 0]) 
and r\ (solid circles, this analysis) in MB Au+Au collisions. 
Error bars include statistical and pr-uncorrelated systematic 
errors, bands show pT-correlated systematic errors. The pair 
of bands at _Raa=1 are the absolute normalization error for 
p+p (larger, dark) and Au+Au (lighter) for it (left) and r\ 
(right). 



The nuclear modification factor Raa is defined as 
l/N cvt dN/dydp T 



Raa 



(Tab) d<jpp/dydp T 



where a pp is the production cross section of the parti- 
cle in p+p collisions, and (Tab) is the nuclear thick- 
ness function averaged over a range of impact parame- 
ters for the given centrality, calculated within a Glauber 
model 16]. When calculating Raa, the measured p+p 
points are used. Raa for r\ production is shown in Fig. [3] 
for four centralities, along with linear fits to Raa- Fit 
parameters are listed in Table HID In the measured pr 
range we observe strong suppression in all but the most 
peripheral collisions. As shown in Fig. 2J for the mini- 
mum bias case the suppression is quite comparable to the 
one observed for tt°, and above 13 GeV/c the (relative) 
systematic errors are smaller. 

TABLE III: Parameters from linear function fit to r\ Raa- 



Centrality 



part 



Slope 



X 2 /NDF 



0-5% 
0-10 % 
10-20 % 
0-20 % 
20-40 % 
40-60 % 
20-60 % 
60-92 % 
MinBias 



351 
326 
236 
280 
142 
61.6 
102 
11.8 
109 



0.008+0.008 
0.011+0.007 

oi 0+ 009 

010+ 007 
u.uiu_ 006 

0.004+0.010 

oio +0 - 018 

0.005+0.011 
0.056j 
0.006+0.007 



,+0.043 
-0.038 



2.77/7 
9.79/7 
11.7/8 
10.8/7 
15.7/8 
4.64/7 
11.7/8 
1.52/6 
10.1/8 



Based upon the most central (0-5%) collisions in [l(| 
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FIG. 5: (Color online) Top: slopes of the linear fits (like the 
ones shown in Fig. along with the fitting errors. Centrality 
is shown in terms of participating nucleons -/V pal .t. Open sym- 
bols are overlapping, solid symbols are non-overlapping cen- 
trality bins (0-10%, 10-20%, 20-40%, 40-60% and 60-92%). 
Also shown: slope of the linear fit to 0-5% tt° data [l(j, shifted 
for better visibility. Bottom: value of -Raa calculated from 
the fit at 5 GeV/c (blue) and 20 GeV/c (red). 



we found that the tt° Raa is consistent with a com- 
pletely flat pt dependence when fitted in the 5 < pr < 
18 GeV/c region, namely the slope of a linear fit was 
m = 0.0017±SS c/GeV. Fitting the current 77 Raa 
data with straight lines gives the slopes and uncertain- 
ties listed in Table Hill and shown in Fig. [5] where central- 
ity is expressed in terms of participating nucleons N part . 
All slopes are consistent with zero; the largest deviation 
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FIG. 6: (Color online) One and two standard deviation x 2 
contours of the linear fits to -Raa in Au+Au collisions for 
0-5%, 0-20% and 20-60% centralities. 
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is less than 2a (for the 0-20% centrality bin). One and 
two standard deviation % 2 contours for selected centrality 
bins are shown in Fig.El For 0-5% centrality we repeated 
the linear fits using only the first 3,4, (n — 1) points 
and found that the slope already stabilizes around its fi- 
nal value with the first few points; data above 10 GeV/c 
improve the significance but barely change the central 
value itself. The same is true for other centralities. 

While the above result indicates that Raa for -q is con- 
sistent with a pr-independent, constant value, and disfa- 
vors a decreasing Raa, a slow rise (~0.01 c/GeV) of Raa 
with increasing pt cannot be excluded. In fact, a detailed 
statistical analysis, comparing to various theories like the 
study done for 7T° in M is necessary once theoretical 
calculations of r\ production are available. However, as- 
suming the linear dependence we can calculate the Raa 
values at 5 GeV/c (where the suppression is already at 
its maximum) and 20 GeV/c; the results are shown in 
the bottom panel of Fig. [5] 

In summary, we measured invariant yields of 77 in 
\J snn = 200 GeV Au+Au collisions at various centrali- 



ties, as well as the rj production cross section in y/s = 200 
GeV p+p collisions in the 5 <pr< 22 GeV/c trans- 
verse momentum range using the PbSc calorimeter of the 
PHENIX experiment at RHIC. The nuclear modification 
factor for 77 in minimum bias collisions is consistent with 
earlier 7r° results. In conclusion, linear fits to Raa as a 
function of pr indicate that Raa is consistent with con- 
stant at all centralities, although a slow rise cannot be 
excluded. 
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